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Abstract 

We use a combination of conventional density functional theory (DFT) and post-DFT methods, 
including the local density approximation plus Hubbard U (LDA+U), PBEO, and self-consistent 
GW to study the electronic properties of Ni-substituted PbTi03 (Ni-PTO) solid solutions. We find 
^ that LDA calculations yield unreasonable band structures, especially for Ni-PTO solid solutions 



S3 
O 

o 



that contain an uninterrupted Ni02 layer. Accurate treatment of localized states in transition- metal 



oxides like Ni-PTO requires post-DFT methods. B-site Ni/Ti cation ordering is also investigated. 
The -B-site cation arrangement alters the bonding between Ni and O, and therefore strongly affects 
^ the band gap (Eg) of Ni-PTO. We predict that Ni-PTO solid solutions should have a direct band gap 

in the visible light energy range, with polarization similar to the parent PbTi03. This combination 



of properties make Ni-PTO solid solutions promising candidate materials for solar energy conversion 
devices. 

PACS numbers: 77.84.Lf, 71.15.-m, 71.20.-b 
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Introduction 

The bulk photovoltaic effect (BPVE) has been observed in ferroelectric perovskite ox- 
ides such as Pb(Zri_ a; Ti :c )03 [1-3], LiNb0 3 [4] and BaTi0 3 [5]. This effect has attracted 
increasing research interest, in view of its potential application for solar energy conversion de- 
vices. As incident light is absorbed by a single-phase ferroelectric material, the photo-excited 
electrons and holes are spontaneously separated, preventing recombination, and producing a 
photovoltaic effect. One ferroelectric material extensively studied for this purpose is BiFe0 3 . 
It has a strong polarization p=* 0.9 C/m 2 [6, 7] and a direct band gap 2.67 eV [8]. Recent 
experiments have demonstrated the BPVE of BiFe0 3 under visible incident light [9-11]. 
However, BiFe0 3 is one of very few ferroelectric materials that have both large spontaneous 
polarization and small E g . Most solid oxide ferroelectrics are wide-gap insulators (E g > 3.0 
eV) that absorb very little of the visible spectrum. 

Recent theoretical work on Pb(M x Tii_ x )03_ x perovskite solid solutions (M=Ni, Pd and 
Pt) provides a promising strategy for designing highly polar semiconducting oxides [12, 13]. 
By substituting the £-site Ti with an O-vacancy-stabilized d 8 M 2+ cation, the resulting 
system displays a decreased E g when compared to PbTi0 3 (PTO), while retaining a large 
polarization. Of these proposed materials, Ni-PTO solid solutions possess great potential 
for practical solar applications, as they can be synthesized from relatively inexpensive raw 
materials. However, due to the strongly correlated Ni-3<i electrons, Ni-PTO has a complex 
electronic structure that is not captured by standard LDA calculations [14]. In this study, 
we explore the electronic structure of Ni-PTO solid solutions with post-DFT calculations, 
including LDA+C/, PBEO hybrid functional and self-consistent GW. We predict the band 
gap and other electronic properties for Ni-PTO solid solutions of various compositions and 
-B-site orderings, to guide future experimental synthesis and measurement. 

Methodology 

We perform first-principles calculations with a plane- wave basis set, as implemented in 
the Quantum-Espresso [15] and ABINIT [16] codes. The LDA exchange-correlation func- 
tional is used for structural relaxations, with a 6 x 6 x 6 Monkhorst-Pack k-point grid [17] 
and a 50 Ry plane-wave cutoff. All atoms are represented by norm-conserving [18] opti- 
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mized nonlocal [19] pseudopotentials, generated with the OPIUM code [20]. The electronic 
contribution to the polarization is calculated following the Berry's phase formalism. [21, 22] 
Several supercell configurations are used, including 1x1x2 and 1x1x3 layered structures 
(with .B-cation alternation along (001) direction); rocksalt (Pb2NiTiOs), \/2 x \/2x2 and 
2x2x2 structures (N1-O5 cages are separated by Ti-Os cages). By studying these configura- 
tions, the effects of different compositions and B-site cation orderings can be investigated. 
Post-DFT methods are then applied to these Ni-PTO systems. In the next section, we 
briefly describe each post-DFT method, to facilitate comparison of results obtained with 
each method. 

Results and Discussion 
Post-DFT methods 

Before applying post-DFT methods to Ni-PTO, we study its parent materials, tetragonal 
PbTi03 and cubic NiO, as reference systems to evaluate the performance of post-DFT 
methods for related materials. We do so because it is known that the Kohn-Sham (KS) 
band gap determined by the conventional LDA method is systematically underestimated by 
50-100 % [23, 24]. To improve the quality of the electronic structure description of Ni-PTO, 
especially the treatment of correlation for the localized Ni-3<i orbitals, we employ LDA+f/, 
the PBE0 hybrid functional, and self-consistent GW (sc-GW) methods. 

LDA+U 

The LDA+£7 method is a computationally inexpensive way to improve the DFT-LDA 
band gap. A simplified version of the rotationally invariant formulation of the LDA+U 
method [14] is employed in the present work, where U is assigned to be the value of the 
spurious curvature of the total energy of the system with respect to the variation in the 
number of electrons (rij) in a localized orbital. [25] The effective interaction parameter U 
is determined by the change of orbital occupations rii with respect to an external potential 
aij projected onto the corresponding localized level: 
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FIG. 1: (Color online) The occupation number (m) of localized d-orbitals for Pb and Ti in PbTiOs, 
with respect to potential shift a. Shown as (black) squares are the values obtained using self- 
consistent (scf) relaxation, and shown as (red) circles are those without (1 st ), for constructing x 
and xo respectively. 



U=(Xo 1 -X- 1 U X« = ^T, (!) 

Here the Xij are the elements of the linear response matrix, rij is the occupation number 
of the localized levels at site i and ex, represents the potential shift applied on the localized 
orbital at site j. \ (the screened response matrix) includes all screening effects from the 
crystal environment which is associated with the localized electrons, and xo (the unscreened 
response matrix) contains non-screening effects of the total energy of the non-interacting 
Kohn-Sham states associated with the system. The latter term can be computed from the 
first iteration in the self-consistent cycle. 

Tetragonal PbTiOs is selected as an example for detailed study. A 2x2x2 PbTiOs 
supercell is used, because it is large enough to give a converged linear response calculation 
of U. Fig. 1 displays the calculated U values for Pb and Ti atoms of PbTiOs. In the 
rotationally invariant LDA+t/ method, on-site interaction energy Eu<xU-[\(l-X)]. For fully 
occupied orbitals such as Pb 5d, X— 1, so on-site correction plays no effect. As displayed in 
Fig. 1 (a), the number of electrons on the Pb 5d orbitals will not change with respect to the 
external potential a, so U=0 eV. Ti has partially filled 3d orbitals, so it has a non-zero U 
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TABLE I: Hubbard U terms (in eV) for Ni and Ti in PbTi0 3 , NiO and Ni-PTO solid solutions, 
obtained with the linear-response approach. 
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FIG. 2: (Color online) Valence band maximum (VBM) and conduction band minimum (CBM) 
versus PBE0 functional exact-exchange fraction a for PbTi0 3 and NiO. Vertical lines represent 
the experimental band gaps, illustrating how ao is chosen (see Table II). 

value. 

Following the procedure above, using the self-consistent linear response method, as im- 
plemented in Quantum-Espresso, the Hubbard U terms for Ti and Ni in NiO, PTO, and all 
Ni-PTO solid solutions can be determined (see in Table I). 

PBE0 hybrid density functional with exact- exchange 



DFT-LDA does not include exact exchange. This is unlike the PBE0 functional, which 
includes a linear combination of exact-exchange and GGA-PBE exchange [26]. The electron 
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TABLE II: Band gaps (in eV) of PbTi0 3 and NiO, calculated using LDA, PBEO (with the optimal 
fraction of exact exchange ccq) and sc-GW methods. 



LDA 


PBEO (a ) 


sc-GW 


Exp. 


PbTi0 3 1.46 


3.50 (0.25) 


3.82 


3.60 [30] 


NiO 1.16 


3.98 (0.125) 


4.14 


4.0-4.3 [31, 32] 



correlation term is the same as that of PBE functional [27]: 

E™ m = <*E X + (1 - a)E™ + E™. (2) 

where a is the fraction of exact exchange. As PBEO includes a portion exact exchange, 
it can greatly reduce the self-interaction error of LSDA or GGA functional, and therefore 
greatly improves the calculated band gap for most semiconductors and insulators. [27] 
However, due to the long-range nature of the exact-exchange interaction and resultant large 
computational cost, (about one to two orders of magnitude more than LDA in the present 
work) PBEO calculations are mainly applied to periodic systems with small unit cells, such 
as those described in the present study. [28] 

In the original PBEO functional, a=0.25. However, it is argued that there is no universal 
a that is applicable for all materials [29]. In order to predict band gap values for new 
materials more accurately, it is considered useful to select an a that yields accurate band 
gaps for the relevant prototype systems. Fig. 2 shows the a dependence of the conduction 
and valence band edges of PbTi0 3 and NiO. Comparing with experimental band gap values, 
we obtain an optimal a =0.25 for PbTi0 3 and 0.125 for NiO (Fig. 2 and Table II). 

Self- consistent GW 

It is widely believed that the self-consistent GW method is the most accurate way to 
obtain electronic band structures, since it includes both the non-locality and dynamical 
correlations absent in DFT-LDA. Therefore, obtaining the self energy, E, is the key quantity 
of any GW calculation. It was shown by Hedin et al. that the self-energy can be represented 
as a product of the Green's function G and the screened Coulomb interaction W , E = 
iGW. [33] 

6 



Since our system has a complicated electronic band structure, in which the band gap is 
dictated by the bonding interactions between the correlated Ni 3d-states and O 2p-states, 
the self-consistent GW (sc-GW) method probably gives the most accurate quasi-particle 
(QP) electronic structures of Ni-PTO. 

Following the methodology developed by Lebegue et al, [34] we employ the contour 
deformation method [35, 36], as implemented in ABINIT, to perform the numerical integral 
for the self-energy. Within this scheme, the contour of the frequency integral for self-energy 
is deformed into integrals along the imaginary and real axes. The numerical integral along 
the imaginary axis can be evaluated by Gaussian quadrature. We find that integrating up to 
114 eV and using 10 quadrature points yields converged results for the imaginary axis. The 
real axis integral is calculated by summing values of the Coulomb screening at a uniform 
mesh of frequencies, with a mesh spacing of 1.25 eV, from to 120 eV. 

This upper limit of real axis integration is chosen from spectral function data, [35, 36] 
(electron energy loss data in the present work), since the correlation part of the self-energy 
can be represented by a spectral function. 

Due to the expensive computational cost (three to four orders of magnitude more than 
DFT-LDA in the present work), we restrict sc-GW calculations to only one Ni-PTO system 
- rocksalt structure. For this configuration, 400 bands are included, with a plane-wave cutoff 
of 25 Ha for the QP and LDA states, and a 2 x 2 x 2 Monkhorst-Pack k-point grid. [17] 
As LDA calculations can qualitatively predict a correct semiconducting ground state for 
rocksalt Ni-PTO (Fig. 6), we choose to represent the initial QP wave functions in the basis 
set of LDA wave functions. QP band gaps are converged to within 0.05 eV. 

Comparison of results for parent end members 

Fig. 3 shows the projected density of states (PDOS) for tetragonal PbTi0 3 and cubic 
NiO, calculated with the LDA, PBE0 and LDA+£7 methods, from which we can evaluate the 
accuracy of each post-DFT method. Experimental studies show that tetragonal PbTi0 3 is 
an insulator with an indirect band gap (X— >T) of 3.60 eV. [30] LDA severely underestimates 
this band gap value by more than 50%. The LDA+£7 method gives a better band gap value 
than LDA, but is still lower than PBE0, which comes closest to experiment. This is because 
the exact-exchange in PBE0 is applied to all states, whereas the Hubbard U of LDA+£7 
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FIG. 3: (Color online) Band gap values and orbital-resolved PDOS of tetragonal PbTi03 and cubic 
NiO using LDA, PBEO and LDA+C7. For PBEO calculations, the optimal ao values are used. 

acts directly only on the localized Ti-3d orbitals. Also, when compared to PBEO, LDA+[7 
significantly alters the PDOS for PbTiC>3, especially for Ti-3<i and 0-2p states. 

NiO is an anti-ferromagnetic insulator with an experimentally measured band gap be- 
tween 4.0 and 4.3 eV [31, 32]. As shown in Fig. 3, LDA predicts a band gap of 1.16 eV and 
LDA+[7 predicts 3.07 eV, while PBEO (with ao=0.125) reproduces the experimental band 
gap. For both LDA and PBEO functionals, the PDOS has a clear split between majority spin 
Ni t2 9 and e g orbitals, while the band gap of the system falls between occupied minority-spin 
t2 g and empty e g orbitals. Applying the Hubbard U to the Ni-3d orbitals increases the NiO 
band gap, relative to LDA. However, LDA+[/ predicts Ni-d spectral weight located at a 
relatively low energy (2-8 eV below Ep). This is in contrast to both experimental pho- 
toemission spectra [37] and theoretical dynamical mean field theory (DMFT) results [38]. 
These studies suggest that there should be a strong Ni-d contribution to the valence band 
of the system. Though LDA+U does predict a larger band gap than DFT-LDA, it does not 
reproduce the salient aspects of the experimental valence-state spectrum of NiO. 

Electronic properties of Pb2NiTiOs solid solutions 

We begin our study of Ni-PTO with the smallest supercell configuration: Pb2NiTiOs. 
In this solid solution, half of the Ti 4+ cations are replaced by O-vacancy-stabilized Ni 2+ 
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cations [12]. This polar substitution will generate a V Q ' and Ni^ ; defect pair (in Kroger- 
Vink notation [39]), where the lowest-energy configuration corresponds to the local dipole 

— * 

parallel to the overall polarization (P). The Ni-PTO structure is tetragonal, with a large 
polarization similar to the parent PbTi03 (Table IV). After replacing Ti 4+ with Ni 2+ , the 
remaining apical O will move away from Ni 2+ , making Ni coordination nearly square planar. 
Since this O atom shifts opposite to the cations, it further increases P. For Pb 2 NiTi0 5 , two 
different B-site cation orderings are studied: (i) 1x1x2 layered and (ii) rocksalt Ni-PTO. 

Layered Ni-PTO is structurally similar to the planar nickelate SrNi0 2 . A previous the- 
oretical study of SrNi0 2 predicted it to be a diamagnetic insulator with fully occupied d xz , 
d yz , d xy and d z 2 orbitals, leaving d x i_ y i empty. [40] Accordingly, since layered Ni-PTO is 
a nickelate with planes of Ni 2+ ions, we anticipate that it should have an insulating ground 
state. 

Fig. 4 (a) shows the band structure for layered Ni-PTO along high-symmetry directions in 
the Brillouin zone (BZ). The LDA band structure reveals a two-dimensional character: the 
energy bands are highly dispersive in the k x -k y plane, but show a very weak z-axis dispersion 
along T-Z. One strongly dispersive band (dashed line in Fig. 4 (a)) crosses Ep, making the 
system metallic in LDA. This band is chiefly composed of Ni-3d x 2_ y 2 and 0-2p. Isosurfaces 
for this band at T, X and M points are displayed in Fig. 4 (b). There is a pronounced 
hybridization between Wi-d x 2_ y 2 and O-p, as the lobes from Ni-d x 2_ y 2 are directed towards 
the four O ions. Specifically, the wave functions at X and M consist of 0(p x )-Wi(d x 2_ y 2)- 
0(s+p z ) and 0(p x )-Ni(d x 2_ y 2)-0(p y ) orbital combinations respectively, both of which display 
anti-bonding character. However, there is a strong bonding profile between Wi-d x 2_ y 2 and 
0-p z at the T point (Ni red (blue) to O red (blue) lobes), which indicates that electron 
density is shared between 0-2p and Ni-d x 2_ y 2 orbitals, again according to LDA calculations. 

LDA calculations show that for layered Ni-PTO, the band gap is zero. We find that in 
LDA calculations, the energy level of the Ni-d x 2_ y 2 orbital is quite low, permitting this low- 
lying orbital to backbond with 0-p z orbitals. This enhancement of the bonding interaction 
between Ni-3d x 2_ y 2 and 0-2p causes LDA to predict a metallic ground state for layered 
Ni-PTO. 

In our LDA+f/ calculations, applying the Hubbard U term (determined to the Ni-3<i 
orbitals, in a self-consistent linear response calculation [25]), yields a split between occupied 
and unoccupied Ni-3d states. This provides a better description of E g for Ni-PTO, relative to 
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FIG. 4: (Color online) (a) Atomic structure of layered Pl^NiTiOs, and its DFT-LDA calculated 
electronic band structure. The band that crosses the Fermi level is shown as a dashed line, (b) 
The wave function isosurfaces are plotted for this band at T, X and M. For clarity, only the 
isosurface near the N1-O2 plane of the band indicated by the dashed lines in (a) is plotted. The 
isosurface values are ± 0.04 e/A 3 (red and blue surfaces). The plots illustrate which atomic orbitals 
compose each eigenstate: T: 0(p z )-Ni(d x 2_ y 2)-0(p z ) , X: 0(p x )-Ni(d x 2_ y 2)-0(s+p z ) and M: 0(p x )- 
Ni(d x 2_y2)-0(p y ) . This series of images shows the drastic change in bonding from T to X to M. 
Pb, Ti, Ni and O atoms are represented as gray, green, cyan and yellow spheres respectively. 

LDA. However, the PBEO hybrid density functional, which includes exact-exchange energy, 
can significantly improve the description of E g of ferroelectric oxides [41] . To obtain the most 
realistic results from a PBEO calculation, the optimal exact-exchange mixing fraction a from 
relevant prototype systems should be used. Unlike PbTiOs, with E g between occupied 0-2p 
and empty Ti-3<i states, both Ni-PTO and NiO have E g between occupied and unoccupied 
Ni-<i states. Therefore, we use NiO as our prototype parent system of Ni-PTO, and select 
a=0.125 for PBEO calculations of Ni-PTO. A comparison of the eigenvalues at T obtained for 
the Ni 3(i-states using LDA, LDA+U and PBEO is presented in the supplementary material 
(EPAPS number). 

Though the most computationally demanding of our post-DFT methods, the many-body 
GW approach is considered the most reliable first principles method to obtain electronic 
band structure, and it has been successfully applied to many kinds of oxides. [35]. In our 
present work, we also employ this method for the ground state calculation of rocksalt Ni- 
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FIG. 5: (Color online) PDOS for Ni-3d and corresponding wave function isosurface for Ni-d x 2_ y 2 
orbital (at T) in layered Pt^NiTiOs, calculated using LDA, LDA+t/ and PBEO hybrid density 
functionals respectively. Compared to LDA, LDA+C/ and PBEO predict smaller downward red 
(blue)-lobes around planar O atoms (see arrows), indicating there is more O s-p hybridization and 
a weaker bonding interaction between ~Ni-d x 2_ y 2 and 0-2p. 

PTO, obtaining a band gap of 1.83 eV. PBEO predicts a value of 1.71 eV, improving upon 
the results obtained by LDA (0.69 eV) and LDA+U (1.59 eV), much closer to the sc-GW 
result. 

Fig. 5 (a) shows the PDOS of Ni-3d in layered Ni-PTO, from LDA, LDA+U and PBEO 
calculations. In both the LDA+U and PBEO PDOS, the Ni-d x 2_ y 2 band lies above Ep, 
leading to the weaker bonding between Ni-d x 2_ y 2 and 0-2p (Fig. 5(b)). Unlike LDA, which 
yields a metallic ground state for layered Ni-PTO, LDA+U and PBEO predict an intrinsic 
semiconductor. Since the LDA+?7 and PBEO results agree and are more consistent with 
the Ni-d orbital paradigm in diamagnetic insulating SrNi02, [40] we consider LDA+U and 
PBEO to be more accurate than LDA, and predict that layered Ni-PTO is a semiconductor. 

The calculated PDOS shown in Fig. 5 also shows other prominent features. As a result 
of the hybridization between Ni and its planar coordinated O, the occupied Ni-c? orbitals, 
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such as d xz , d yz and d z 2, have a large band width of 2 eV (both LDA and PBEO results). As 
the Ni ion is shifted up, away from the oxygen plane, the hybridization between Ni-rf z 2 and 
0-p z can give rise to bonding and antibonding states. The double peaks in d xz , d yz states 
should correspond to bonding and antibonding states, which come from the hybridization 
between Ni-d xz (d yz ) and 0-p x (p y ) (for details see the supplementary material). The strong 
hybridization between Ni and planar O ions indicate the covalency of Ni-0 bonds in layered 
Ni-PTO. 

After exploring the electronic properties of the layered structure, we extend our investi- 
gation to rocksalt ordered Ni-PTO (Fig. 6 (a)). In the layered structure, there is a strong 
bonding character between Ni-d x 2_ y 2 and 0-p z (at T). However, in rocksalt Ni-PTO, there is 
an anti-bonding relation between Ni-d x 2_ y 2 and 0-p z , while bonding with 0-p z occurs almost 
exclusively between Ti and O (Fig. 6 (b)). This hybridization pattern indicates that the 
electron density is mainly shared between Ti-3d and 0-2p, and that the Wi-d x 2_ y 2 orbital 
is almost empty. This leads to the prediction of a semiconducting ground state for rock- 
salt ordered Ni-PTO, even at the LDA level (Fig. 6 (c)). A consistent picture emerges, as 
LDA+C/, PBEO and sc-GW all predict rocksalt Pb 2 NiTi0 5 to be a semiconductor. 

Besides the Ni-d x 2_ y 2 orbital, other Ni-rf bands also differ between the rocksalt and layered 
structures of Ni-PTO. In layered Ni-PTO, a strong hybridization between Ni and O yields 
energy bands with a large band-width. However, in rocksalt Ni-PTO, interrupting the Ni0 2 
layer with Ti-0 cages makes it more favorable for the filled 0-2p orbitals to mix with the 
empty 3d orbitals of Ti 4+ [42]. Therefore, Ni-3rf orbital is weakly hybridized with 0-2p, 
reducing its band width in rocksalt Ni-PTO (Fig. 6 (c)). 

i?-site cation ordering effects 

The results of our electronic band structure analysis indicate that different B-site cation 
orderings for Pb 2 NiTi0 5 systems can give rise to different degrees of Ni-0 bond covalency. 
To support this, we calculate the Born effective charges (BECs) for £-site cations of Ni-PTO. 

Since LDA incorrectly predicts a metallic ground state for layered Ni-PTO, the BECs 
and polarization P of Ni-PTO should be calculated at the post-DFT level. The BEC tensor, 
Z* a/3 , is defined as the change in the macroscopic polarization along direction a in response 
to a rigid displacement of a sublattice of atoms k in direction (3, times the unit cell volume 
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FIG. 6: (Color online) (a) Atomic structure of rocksalt -B-site ordered Pl^NiTiOs. A 2x2x2 
supercell containing 36 atoms is shown for clarity, (b) Isosurface of LUMO band (LDA results) at 
r. It is mainly composed of d x 2_ y 2 orbitals, displaying anti-bonding character between ~Ni-3d x 2_ y 2 
and 0-2p z . (c) PDOS for Ni-3d orbital, obtained using LDA, LDA+C/ and PBEO functional 
calculations respectively. 



n [43]: 



dP a AP a 

Z *K,a/3 = ^0^— |^ =0 - fi 7T— \e=0 ( 3 ) 



K,fi 



Following this definition, Z* can be estimated by the finite-difference method: each atom 
k is displaced slightly from its equilibrium position by Ar^, and the resulting change in 
polarization AP a can be calculated using the Berry-phase method. [21, 22] The displace- 
ment should be small enough to ensure the validity of the linear treatment of Eq. (3). For 
each structure, the U value is determined by a linear-response calculation, so that LDA+U 
calculations can generate an insulating ground state for all Ni-PTO structures. In this way, 
we calculate both Z* and P at the LDA+U level. 

Table III gives the calculated non-zero BECs for Ni ions in Ni-PTO systems. As Ni 2+ ion 
strongly favors a four-fold coordinate environment, its planar BEC tensor elements {Z* x and 
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TABLE III: The Born effective charges Z* of Ni cations in Ni-PTO solid solutions. Z* is calculated 
using LDA+C/ finite difference methods. 



Composition 


*v* y* 

yy zz 


1x1x2 Pb 2 NiTi0 5 
1x1x3 Pb3NiTi 2 08 


2.29 2.29 0.44 
2.33 2.33 0.37 


rocksalt Pb 2 NiTi0 5 
v / 2x v / 2x2 Pb 4 NiTi 3 0n 
2x2x2Pb 8 NiTi 7 2 3 


2.14 2.14 0.86 
2.17 2.17 0.72 

2.15 2.15 0.78 



Z* y ) are larger than its nominal charge, indicative of covalency for the planar Ni-0 bonds. 
Since the apical O moves away from Ni to create a square planar environment, the bonding 
between Ni and the apical O is weak. In the finite-difference simulation of Z*, when Ni ion 
moves closer to the apical O, the electron density will transfer from Ni to O. Because of this 
negative charge density flow, the Ni ion will have a Z* z smaller than its nominal charge. 
Moreover, the following trend is identified: as .B-site cations change from layered to rocksalt 
ordering, there is a decrease of the planar BEC tensors of Ni, and an increase of Z* along 
Ni-0 vacancy direction. This trend is in good agreement with the bond covalency results 
we obtained above, that layered Ni-PTO has a more covalent Ni-0 bond and larger Z* x and 
Z* y tensor for Ni, while less covalent Ni-0 bond in rocksalt structure leads to the smaller 
in-plane Ni BECs. 

In perovskite AB0 3 compounds, a more ionic B-0 bond usually yields a larger E g , while 
a more covalent B-0 bond is associated with a smaller E g . This trend has been verified for 
many perovskite systems, such as Pb(Zr :E Tii_ :E )03 [44] and sulfide perovskite BaZrS 3 [45]. 
In light of this trend, we now examine the Ni-PTO solid solution systems. 

E g of several Ni-PTO solid solutions are summarized in Table IV. Relative to PTO, E g of 
Ni-PTO solid solutions are considerably reduced. For the rocksalt structure, LDA+£/, PBE0 
and sc-GW calculations increase E g (relative to LDA) by different amounts, and the trend 
in calculated E g is: LDA<LDA+t/<PBE0<sc-GW. Unlike the parameter-free sc-GW, 
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TABLE IV: Reported here are the c/a ratio, polarization P (in C/m 2 ), as well as indirect Eg (in 
eV) calculated using different methods for all the Ni-PTO solid solutions under investigation. Also 
shown are tetragonal PbTiC>3 results. For Ni-PTO solid solutions, direct band gaps at T are given 
in parentheses. Structural parameters are obtained by LDA calculations and, P is calculated using 
the LDA+C7 finite difference method. We use the optimal a for PBEO calculation. 



Composition 


c/a 


P 


"kg 


E LBA+U 


I7.PBE0 
-^g 


E sc- GW 


PbTiOg 


1.05 


0.81 


1.46(3.05) 


2.17(3.80) 


3.50(4.31) 




1x1x2 Ni-PTO 
1x1x3 Ni-PTO 


1.14 
1.14 


0.85 
0.84 


0.00(0.00) 
0.00(0.00) 


0.70(1.19) 
0.52(1.57) 


0.87(1.02) 




rocksalt Ni-PTO 
y/2xy/2x2 Ni-PTO 
2x2x2 Ni-PTO 


1.21 
1.12 
1.08 


0.96 
0.90 
0.83 


0.69(1.12) 
0.84(0.97) 
0.70(0.94) 


1.59(2.65) 
1.53(2.24) 
1.60(1.90) 


1.71(2.64) 


1.83(2.70) 



E g results obtained from LDA+£7 and PBEO fairly strongly depend on the computational 
parameters (exact-exchange mixing fraction a for PBEO and Hubbard U term for LDA+£7). 
In the present work, LDA+[7 and PBEO E g are in good agreement with sc-GW results, 
verifying our parametrization for these calculations. It is worth noting that even though 
LDA+£7 yields E g similar to PBEO and sc-GW, LDA+£7 gives a more localized and perhaps 
less realistic DOS distribution (Fig. 5 (a) and Fig. 2 in supporting information). Therefore, 
we only use it to obtain the relative E g trend, but do not rely on LDA+£7 for the precise 
electronic structure. Our calculated E g exhibits obvious 5-site ordering effects: layered Ni- 
PTO always has a smaller E g than rocksalt ordered Ni-PTO, and the variation in Ni doping 
percentage has little impact on E g value. 

By comparing the BEC tensor and E g results of Ni-PTO solid solutions, we can establish 
the relationship between B-site cation ordering and band gap. The electronic properties of 
Ni-PTO primarily depends on the Ni-0 bonding. In layered Ni-PTO, there is a complete 
Ni-O-Ni network, and strong hybridization between Ni-d and O-p increases Ni-0 bond cova- 
lency, yielding a smaller E g . However, if the Ni-O-Ni network is interrupted by Ti-0 bonds, 
the hybridization between the filled 0-2p orbitals and the empty 3d orbitals of Ti 4+ will be 
more favorable [42]. As a result, there is a more ionic bond between Ni-o? and O-p, and the 
system has a larger E g . 



15 



Guided by the E g trend obtained above, we estimate that the lower E g limit of Ni-PTO 
solid solutions should come from the layered structure, which could potentially be synthe- 
sized via layer-by-layer deposition techniques; the upper E g limit corresponds to the more 
ionic rocksalt ordering structure, which can perhaps be obtained with a thermal annealing 
process. According to our prediction, layered Ni-PTO should have a direct optical band gap 
comparable to bulk Si, and E g of rocksalt ordered Ni-PTO is close to that of BiFe0 3 [8]. 
In this way, Ni-PTO solid solutions should have direct band gaps in the visible light energy 
range, displaying semiconducting electronic properties. 

Conclusion 

In conclusion, we present theoretical studies on the electronic properties of 
Pb(Ni 2 .Ti 1 _ 2 .)0 3 _ x solid solutions. We identify the relation between 5-site cation order- 
ing and E g of Ni-PTO systems. We expect that this relation can be used a guiding principle 
to understand E g behaviors for many perovskite solid solutions. Based on our theoretical cal- 
culations, we predict that the Ni-PTO solid solution can offer both a large polarization and 
a direct band gap in the visible light energy range. This combination of properties indicates 
that high photovoltaic efficiency can be realized in solar devices based on Pb(Ni a ,Ti 1 _. E )03_. E 
oxides. 

Furthermore, we demonstrate that the standard DFT-LDA provides qualitatively incor- 
rect electronic structure for layered Ni-PTO, and potentially other systems being investi- 
gated as possible components in solar harvesting devices. A more accurate treatment of 
the Ni-3cZ states of Ni-PTO required post-DFT LDA calculations to obtain reasonable band 
gaps. Of the work presented here, the LDA+£7 method is the least computationally ex- 
pensive of the post-DFT methods, though it sometimes yields unreasonable crystal field 
splits and unphysical PDOS features. Though the most computationally demanding, sc- 
GW yields band gap results that are closest to experiment. The band gap results of PBEO 
are close to that of sc-GW and reproduce salient spectral features of the PDOS, while the 
computational cost is less than that of GW. For the purposes of calculating accurate band 
gaps, we suggest that the PBEO method is a reasonable compromise of accuracy and cost. 
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Supplementary materials - Post density functional theoretical 
studies of highly polar semiconducting Pb(Tii_ a; Ni a ;)03_ a ; solid 
solutions: The effects of cation arrangement on band gap 
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FIG. 1: (Color online) The wave function isosurfaces of the Ni-3d orbitals at T from the layered 
Pb2NiTi05, calculated by LDA. The isosurface value is ± 0.04 e/A 3 . Except for the non-bonding 
Ni-3d xy state, the wave functions show obvious hybridization between Ni-3d and 0-2p. a and a* 
correspond to the 7r-type bonding and antibonding states between Ni-3d 2 2 and 0-p z . b and b* are 
degenerate bonding and antibonding states from Ni-d xz (d yz ) and 0-p x (p y ). 
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FIG. 2: (Color online) The energy levels (in eV) for the corresponding eigen states in Fig. 1, 
calculated by LDA, PBEO and LDA+J7 methods. Fermi energy level is taken as energy zero. Both 
LDA+C/ and PBEO predict an empty Ni-d a ,2_ y 2 state. However, LDA+C7 predicts energy levels of 
occupied Ni-d states far below the Fermi level, leading to the overall improper PDOS feature (Fig. 
5 (a) in the main text) for layered Ni-PTO. Note that PBEO has similar PDOS to LDA, and has 
d x 2_ v 2 unfilled above Ep. 
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